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ABSTRACT 


The Optical cavity of a Free Electron Laser (FEE) is composed of components that must 
be maintained to very tight tolerances. The shipboard environment is one that will 
preclude a direct coupling of FEE components to the ship. This thesis will explore the 
basis for these tight tolerances, and how to isolate them from the FEE. 

A solid model of a potential FEE system will be developed using SolidWorks. 
This model will then be converted to a finite element model in ANSYS. The finite 
element model will be used to calculate the system’s eigenvectors. These eigenvectors 
will be used to develop a state space model in MATLAB. Driving functions simulating 
sea state 6 and under water explosion will then be applied to the state space model and 
the motion of various components will be tracked. This simulated motion will be used to 
develop and test a passive control system to damp out the vibrational input to the FEE. 

It is not possible for a passive system to totally isolate the FEE from excitation by 
the ship environment. A passive system that minimizes the inputs to an active control 
system will be developed. An active system that will handle the final mirror stabilization 
for a FEE optical cavity will be left for further research. 
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I. INTRODUCTION 


A. BACKGROUND 

There are several distinet eombat seenarios where the free eleetron laser (PEL) 
will be of great benefit to the Navy. For example, the FEE can become a ship’s primary 
self-defense weapon against inbound airborne threats, but can also provide a less-than- 
lethal means of deterring small craft, and can be used for a precision strike against shore 
targets near the coastline. 

1. Free Electron Laser for Self Defense 

a. Aircraft 

The FEE has the capability to project energy long distances at the speed of 
light. A megawatt class FEE will also be capable of destroying enough material in about 
one second to render a target missile incapable of flight. Other types of aircraft will 
require more or less time on target depending on their construction. The laser is also 
capable of retargeting nearly instantly since the beam director need only be aligned to the 
new target. Unlike some other high power lasers, the FEE does not have a required 
cooling time between or after a certain number of shots. 

b. Waterborne Threats 

The optical targeting of the FEE is expected to be accurate enough to 
enable it to be used against small craft in a number of ways. It can be used at lower 
power settings to cause small amounts of damage to the craft as a sort of a warning shot. 
It also can be used to make holes in the hull near the waterline or to disable the engines. 

2, Shore Engagements 

Some targets on shore can be disabled or destroyed by the loss of a small amount 
of material. Since the operation of the FEE only requires electricity, the only cost in its 
use is ship’s fuel. The “pinpoint” accuracy of the system allows for engagements that 
would have been unfeasible with any other system based on a cost benefit analysis. 
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B, SENSITIVE COMPONENTS 

The many benefits of the PEL system eome at a non-trivial cost. The system is 
costly but it is not likely to be any more so than other weapon systems of equal power. 
Any system that works with megawatts of power is composed of many sensitive 
components with tight tolerances; the FEE is no exception. The laser is quite large and 
has components that must be aligned within microradians and must also have their 
separation maintained within micrometers. The operating environment may be one 
which vibrates with significant amplitude over a broad spectrum of frequencies. 
Achieving either of these conditions would be no small feat in a shipboard environment; 
the FEE requires both. 

C. SCOPE OF RESEARCH 

The system to maintain tolerances on the FEE will necessarily be one that is quite 
sophisticated. The scope of this thesis will be limited to the alignment of three important 
components: the two optical resonator mirrors, and the electron beam. This project will 
model the whole FEE system, simulate that model in a shipboard vibrational 
environment, develop a passive isolation system to minimize vibrational input, and 
establish the boundary conditions for an active control system to stabilize these three 
components. 
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II. FREE ELECTRON LASER SYSTEM 


A, COMPONENTS 

The Free Electron Laser (EEL) is comprised of several main types of components. 
These components range from bending magnets to superconducting linear accelerators. 
A complete assembly of these components is in Eigure 1. The figures in this section are 
for illustrative purposes only; this model is a notional design that captures the mechanical 
properties of various EEL components for the purpose of studying the effects of a 
vibrating environment on the system. 



Resonator Mirror 


Quadrupole Set 


Beam Dump 


Accelerator 


Undulator 


U-Turn 


Dipole 


Injector 


Figure 1: FEL Solid Model 


3 





































1, Injector 

The injector, an example of which is in Figure 2, is a source of electrons. In this 
case they are generated by a short pulse length laser that is impinging on a photo cathode. 
The laser pulses used to generate the photoelectrons in the photocathode are synchronized 
with the radio frequency (RF) field in the injector cavity to produce positive acceleration 
to the electrons. These electron bunches are also produced in phase with the RF Field in 
the main accelerator so that their energy can be further increased. The injector requires 
vacuum, RF power, and cooling systems in addition to the seed laser to support its 
operation. The injector, as well as the entire beam line, requires extremely high vacuum; 
approximately 10-**’ Torr is necessary to prevent unwanted electron scattering, and 
cathode degradation. 


three cell 
RF cavity 


Ultra precise 
photo cathode 


Figure 2: Injector Cutaway. From [1] 

2, Accelerator 

The accelerator uses an RF power source to establish an oscillating 
electromagnetic field inside a series of cavities such as those pictured in Figure 3. This 
field will incrementally, depending on phase relation, accelerate or decelerate a bunch of 
electrons as they pass through each cavity. The phase relationship between acceleration 
and deceleration allows for energy recovery with a recirculating beam line. In this 
model, the accelerator cavities are superconducting devices operating at 2K. Working at 
this temperature requires additional support systems using liquid helium. The thermal 
insulation for this device will be accomplished in steps. The 2K volume will be 
surrounded by an insulating vacuum, outside this vacuum layer is a 77K liquid nitrogen 
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layer, next is another insulating vaeuum, finally a typieal thermal insulation is used. The 
insulating vaeuum is not required to be of the same quality as the beam line and is only 
held at approximately 10-"^ Torr. The aeeelerator eavities however, being part of the 
beam, eontain the same 10-'*’ Torr vaeuum as the injeetor. 



Figure 3: Accelerator Cutaway. From [2] 


3, Bending Magnets 

There are several types of bending magnets that may be used in the PEL. In this 
design, there is a smaller dipole merge magnet, and a much larger magnet for turning the 
beam around creating the recirculating loop. 
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a. 


Dipole 



Figure 4: Cross Section View of a Dipole Magnet with Attached Piping and 

Windings Shown in Brown 

The smaller dipole magnets eause a bend radius that is proportional to 
eharge, energy, and magnetie field. A cross section view of this type of magnet can be 
seen in Figure 4. In the figure, the grey material is steel and the brown is copper 
simulating the windings for this electro-magnet. These colors are constant throughout the 
SolidWorks [3] model and will not be described again. These magnets are used in 4 
places in this FEL design: (1) They interlay low and high energy electron bunches at the 
outlet of the injector so that they enter the accelerator along a common beam line axis but 
at different RF phases, (2) they separate low and high energy electron bunches at the 
outlet of the accelerator so that the high energy electrons travel down the main beam line 
and the low energy ones enter the beam dump, they also allow the optical beam and 
electron beam to (3) enter and (4) exit the optical cavity. For this model, the same 
magnets are used with a slightly different piping arrangement to deflect the electron beam 
and, where necessary, to align it for travel through downstream components. 
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U-Turn 


b. 



Figure 5: 


Cross Section View of a U-Turn Magnet with Attached Piping and 
Windings Shown in Brown 


Large 90° dipole bending magnets are used at opposite ends of the eleetron 
beam line to create the loop that allows for energy recovery in a recirculating PEL. A 
cross section view of one of these magnets can be seen in Figure 5. Two of the 90° 
magnets are used in conjunction with a focusing magnet to turn the beam a full 180°. 

4, Focusing Magnets 



Figure 6: Quadrupole Focusing Magnet with Attached Piping and Windings 

Shown in Brown 
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Quadrupole focusing magnets are used to foeus the electron beam radially at 
various loeations around the beam line. These magnets are used in 19 locations along the 
beam line to maintain the approximately 1 millimeter beam envelope radius. This model 
uses these magnets in sets of three to balance the induced angular spreads in both of the 
transverse directions being controlled. An example of one quadrupole magnet can be 
seen in Figure 6. 

5, Undulator 



Figure 7: Undulator Principle. From [4] 


The linear undulator is a device that contains a set of magnets whose poles 
alternate between north and south along the long axis of the undulator. An example is 
sketched in Figure 7. The alternating of poles sets up a sinusoidally varying magnetic 
field that will, via the Lorentz Force, cause the electron beam to travel in a sinusoidal 
path. The acceleration associated with a sinusoidal trajectory will cause the eleetron 
beam to emit synehrotron radiation. 

6. Beam Dump 

The beam dump is a large block of copper surrounded by shielding. Its purpose is 
to capture the electrons that come out of the accelerator at about 5 MeV and then to 
shield the radiation produced by the copper-electron interaction. The beam dump must 
also be eooled as it will need to dissipate about five megawatts of heat generated by 
stopping the electrons. 
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7. 


Resonator Mirrors 


The navy PEL may be an oseillator design, as is considered here. This design will 
use a pair of resonator mirrors set up to reflect light pulses through the undulator in 
synchronism with the electron pulses. Furthermore, one of the mirrors will be partially 
transparent to allow for usable light to “leak” out of one end. These mirrors will need to 
be a significant distance apart to allow for sufficient spreading of the optical beam to 
prevent distorting or damaging the mirrors, ft is this distance and the sensitivity of the 
alignment with the electron beam in the undulator that this project will be exploring. 

B, OPERATION 

Electrons are liberated from the cathode in the injector in bunches with charge of 
about one nanocoulomb at a rate of 750 MHz, resulting in approximately one ampere of 
current. These bunches are then accelerated to about five MeV in the injector, and sent to 
the accelerator. In the accelerator, the electrons enter in phase with the oscillating RE 
field, where their energy is increased to about 100 MeV. The electrons then travel 
through a series of focusing and bending magnets along their path to the undulator. The 
sinusoidally varying magnetic field in the undulator will cause a periodic acceleration due 
to the Eorentz force. The electrons will be directed off on their trajectory in the plane 
normal to the magnetic field in a sinusoidal path; this transverse acceleration causes the 
electrons to interact with the laser field in the resonator and transfer a few percent of their 
energy to the light. The electrons then travel through another series of bending and 
focusing magnets on their trip back to the accelerator. This time, the bunches of 
electrons enter the accelerator 180° out of phase and are decelerated so that the high 
energy electrons give up most of their energy to the RE Field. This both reduces the 
required RE energy input to the accelerator, and drastically lowers the radiation and heat 
produced in the beam dump. 
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Figure 8: FEL Power, Gain, and Optical Spectrum vs. Frequency for 2500 Passes 

through the Undulator. From [5] 


This interaction allows the light to become coherent at a specific frequency over a 
large number of passes through the undulator. The interaction occurring many times also 
allows for high power to be developed in the optical beam. Figure 8 demonstrates the 
development of the optical mode in this type of FEL. The typical gain spectrum is shown 
in purple, as a function of wavelength, with higher gain above the zero line. The gain 
spectrum plot also shows the threshold for continuous operation in yellow, gain above 
this line is necessary to make up for the energy extracted. The blue plot shows the 
evolution of the power spectrum of the optical mode for each pass over 2500 passes with 
darker blue indicating more optical power. Finally, the power spectrum at the end of the 
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2500 passes is shown in black. The simulation is started with the spontaneous optical 
spectrum (red) and the spectrum rapidly narrows to a very small bandwidth. This process 
happens quickly in the time domain: 


However, the electrons must be able to repeat the process with high accuracy and 
repeatability several thousand times. 



Figure 9: FEL Extraction (ti) vs. Desynchronism (d). From [5] 

The light produced in the undulator is stored between the two resonator mirrors, 
and is amplified each pass through the undulator. In order for this amplification to occur, 
the distance between the mirrors must be set so that the arrival of the optical pulse is 
synchronized with the arrival of an electron pulse so that they enter the undulator 
together. Figure 9 shows a plot of FEL extraction percentage (r|) vs. the range of mirror 
positions AS quantified by the dimensionless desynchronism parameter 
One of the mirrors is moved closer to the undulator by the distance AS in order to account 
for optical lethargy [5]. This ensures that the two pulses arrive in synchronism many 
times and allows for coherence and high power to be developed in the optical beam. 
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From Figure 9, the range of possible values for d is about 0.45 in this example. 
Assuming an optieal wavelength of A = 1 p,m, and the number of periods in the 
undulator is N = 20, the mirror ean ehange by about AS = d * NA/2 = 0.45 * 20p,m/2 = 
4. 5p,m before this representative system stops working. 

C. SYSTEM TOLERANCES 

In the shipboard environment, vibrations from various sourees ean eause 
fluetuations in the separation and alignment of the resonator mirrors. From the 
ealeulation of AS above, one of the toleranees on the position keeping of the resonator 
mirrors ean be determined. Moving one mirror through a range of Sfim takes the laser 
from not working, to produeing max power, to not working again. This large a swing in 
output power is not aeeeptable for FEL operation; the extraction needs to be better 
controlled. Because both mirrors move, each must be maintained within about 1/rm to 
maintain approximately constant extraction. Setting the starting position of each mirror 
at the center of its band, the tolerance on position in this direction is ±0.5/rm. 
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Mirror Tilt-JLab FEL 


3 
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0 5 10 

AO^ (ill rad) 


Figure 10: FEL Extraction vs. Mirror Tilt Study Performed on a Power Upgrade 

to the JLab FEL. From [6] 

A study of mirror tilt was also performed at NPS on a Thomas Jefferson National 
Accelerator Facility’s (JLab) FEL [6] and the results are shown in Figure 10. The top 
sketch shows 0m as the angle between the long axis of the undulator and the mirror being 
rotated. The middle plot shows the electron beam (red), the electron pulse shape 
(purple), the optical mode at the beginning and end of the undulator (yellow), and the size 
and shape of the resulting optical beam (light blue) for a deflection of 5 microradians. 
The resulting FEL extraction is shown in the bottom of Figure 10 as a function of 0m. Of 
note, this study showed that a mirror tilt of only 5 microradians produced a tilt in the 
optical mode of 1 milliradian and a reduction of laser power by a factor of two. This 
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serves to emphasize why eontrolling the mirror angle so tightly is both important and 
neeessary for PEL operation. These results establish two new toleranees on the resonator 
mirrors that will be explored in this projeet, the angle of the two resonator mirrors must 
be maintained within about ±5 mieroradians, and the optieal mode rotation must be 
maintained within about ±1 milliradian. 


( H=^ 



Mirror Shift-JLab PEL 



O'-'-'- 

0 50 100 150 

(pm) 


Figure 11: FEL Extraction vs. Mirror Shift Study Performed on a Power Upgrade 

to the JLah FEL. From [6] 


To prediet the effects of mirror shift, another study was performed on JLab’s FEL 

in [6], as shown in Figure 11. The top sketch shows AYm as the distance between the 

long axis of the undulator and the center of the mirror being displaced. The middle plot 

shows the electron beam (red), the electron pulse shape (purple), the optical mode 

14 














(yellow), and the size and shape of the resulting optical beam (light blue) for a deflection 
of 75 micrometers. The resulting PEL extraction is shown in the bottom of Figure 10 as a 
function of AYm. This 75 micrometer shift also produces an optical mode rotation of 1 
milliradian and lowers extraction by a factor of two. From these results, a tolerance of 
±50 micrometers is established on resonator mirror shift and the ±1 milliradian tolerance 
on optical mode rotation is reinforced. 



Gain vs. electron beam tilt 


50 


O 20 



lasing threshold 


0.5 , J 

0y (mrad) 


1.5 


Figure 12: Gain vs. Electron Beam Tilt Setup. From [5] 


Another similar study of was performed [5] focusing on electron beam tilt in a 
typical FEE as shown in Figure 12. The top sketch shows 0y as the angle between the 
optical mode and the electron beam. The middle plot shows the electron beam (red), the 
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ideal Gaussian shaped optical mode (yellow), the size and shape of the optical beam 
(light blue), for an angle 0y between the optical mode and electron beam. This study 
worked on the premise that prior to entering the undulator, the electron beam was 
deflected out of alignment with the optical resonator. The result was a loss of optical 
weak field gain, shown in the bottom figure, since the two beams were not always 
collocated. If there is no light present near the electron beam, stimulated emission cannot 
occur and there will be reduced gain in that region. Without adequate gain, the PEL will 
be unable to start up. If it is assumed that the electron beam travels through the center of 
the undulator, as designed, then a logical extension of this study is that if the optical beam 
is out of alignment by the same angle as the electron beam in this study, the same results 
will occur. At about one milliradian of optical beam misalignment the laser loses 75 
percent of its weak field gain and becomes unlikely to start up. The mirror alignment 
required for FEE startup is ±0.5 milliradians. 

In [7], a simulation of JEab’s FEE oscillator was performed to determine 
constraints on energy extraction (p) from the electron beam due to its shift (Ay) and tilt 
(AOy) together relative to the undulator axis as shown in Figure 13. Maintaining 
extraction at approximately two percent is estimated as necessary to generate the useful 
optical power that the FEE is intended to produce. For the same reasons as before, both 
the angular and shift results for the electron beam can be applied to the optical beam. 
The optical beam angle must be controlled to within 1 milliradian. Another new 
tolerance on mirror motion emerges from this data, the optical beam and thus the optical 
mode shift must be controlled to ±0.6 millimeters. 

The tolerances developed in this chapter are summarized in Table 1. 
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Figure 13: Extraction vs. Electron Beam Shift and Tilt. From [7] 


Table 1: FEE System Tolerances 


Dimension to be Controlled 

Tolerance 

Mirror Angle 

±0.5 microradians 

Mirror Shift 

±50 micrometers 

Optical Path Angle 

±1 milliradian 

Optical Path Shift 

±0.6 millimeters 

Mirror Separation 

±0.5 micrometers 
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III. MODELING 


A. SOLID MODEL 

The modeling began with the development of a solid model of a notional PEL 
oscillator system using SolidWorks. A tour of the now disassembled Superconducting 
Accelerator (SCA) FEE at Stanford University was taken to familiarize the designers 
with the components of the system, and their construction. 

The first step in the design process was to develop a table of masses of 
components from several existing models; see Table 2. The three columns from 
Advanced Energy Systems (AES) are different power levels that were presented in [5]. 
The FEESIM [8] column is based on a megawatt class laser simulation. Where 
information was not available, masses were developed using common engineering 
materials and approximate dimensions of components. Using the data compiled, and 
standard engineering materials, a solid model of a notional FEE system was developed. 
There are two places where the constructed model for this project differs greatly from the 
data collected. The injector modeled is a superconducting quarter wave gun being 
developed by Niowave, Inc. in Eansing, MI. The second deviation is that this model is 
using a half meter permanent magnet undulator. Both of these designs are much less 
massive than their traditional counterparts. The exact properties of the system with 
regards to the functionality of the laser itself are not critical to this project. What is 
important, and what is being emphasized in this model, is that individual mass, moment 
of inertia, and stiffness for various components be consistent with a possible 
configuration for a FEE system. A three dimensional rendering of the solid model can be 
seen in Figure 1 above. 
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Table 2: 


Component Mass Comparison from Other Models 


Component 

AES 

Model #1 
mass (kg 
total) 

AES 

Model #2 
mass (kg 
total) 

AES 

Model #3 
mass (kg 
total) 

EEESIM 
90kW 
(kg each) 

Model 

Average 

Selected 
Project 
Model 
mass (kg 
each) 

Photocathode 

Laser 

250 

450 

704 

* 

468 

★ 

Injector 

2700 

2900 

4000 

1171 

2693 

173 

Dipole 

* 


* 

37 

37 

50 

Quadrupole 

magnet 

* 

* 

* 

262 

262 

269 

LINAC 

2350 

1900 

1900 

2121 

2068 

1717 

Beam dump 

400 

1300 

2300 

* 

1333 

787 

90° Dipoles 

* 


* 

* 


360 

Undulator 

1000 

800 

800 

317 

729 

100 

Mirror 

Assembly 

=1= 


* 

* 


118 

optical path 

* 

* 

* 

500 

500 

843 

Cryo system 

14000 

13000 

10000 

* 

12333 

★ 

RF Sources 

1100 

8100 

20000 

* 

9733 

★ 

RF 

distribution 

2400 

5400 

13000 

* 

6933 

★ 

HVDC 
power supply 

2000 

4000 

12000 

* 

6000 

★ 

Total System 
Mass (Tons) 

55 

67 

85 

36.2 

61 

40.5 

* Data Not Available * Component not Included in Model 
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B. 


ANSYS MODEL 


In order to facilitate simulations, a finite element model based on the solid model 
needs to be created. A solid model of a system such as this contains a huge amount of 
data. There are techniques to generate a finite element mesh directly from the solid 
model; however, the matrices associated with them become quite large and are prone to 
errors and long processing times. To avoid the long simulation times and possibly 
inaccurate results, a beam element based finite element model was developed using 
ANSYS [9]. The finite element model is one that is effectively composed of sticks and 
masses. This model uses the spatial properties of the solid model and joins them with 
linear elements that have the same mass and inertial properties as the bodies developed in 
the solid model. Top views of both models can be seen below in Figure 14 to illustrate 
both the geometric similarity between the two models and how drastically simplified the 
finite element model (top) is as compared to the solid model (bottom). As drastic as this 
simplification is, there is no loss in the how the model behaves globally, only in how it is 
processed in the simulation. 



Figure 14: Finite Element model (top). Solid Model (bottom). The Details of the Solid Model 

are Shown in Figure 1 
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1, Element Development 

The reason that there is no loss in fidelity with the finite element model is because 
of the way the elements are programmed, and the questions asked of the simulation. 
Each element has two sets of properties associated with it, one for material properties, 
and one for inertial properties. The set of material properties contains three values: 
Young’s modulus, Poisson’s ratio, and the density. The set of inertial properties for the 
structural beam element used in this model contains 12 values: cross-sectional area, area 
moment of inertia in all three directions, thickness in both cross-sectional directions, 
orientation angle of the element, shear modulus, and others that are not used in this 
model. 

These values are obtained in one of two ways. If the component being modeled is 
made of a single uniform material then the general material property values can be looked 
up in a handbook and the other properties can be measured off of the solid model. If the 
component is a composite or has a non-uniform cross section, then an average value must 
be determined for that component. A dimensioned cross section of a quadrupole magnet 
is shown in Figure 15; example composite body calculations of area moments of inertia 
and density follow, using standard formulas and the parallel axis theorem. All 
dimensions in Figure 15 are in meters and directly translate into the formulas following 
the figure. 
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All of the other properties for these elements are similarly calculated. 


2. Eigenvectors and Eigenvalues 

Based on the elements data, and geometry entered, ANSYS produces an output 
file that contains the eigenvalues, and eigenvectors of oscillation for the first 50 
vibrational modes of the system. These eigenvalues are the natural frequencies for each 
of the modes and the eigenvectors are the shapes of each of the modes. This file will then 
be imported into MATLAB for further processing and application of a driving function. 

C. SEA SIMULATION 

The sea input to the model was approached in two ways; a steady state input 
representing a series of waves and a shock input. For the purposes of this project, a 
DDG-51 class ship will be used as the vessel to carry the FEL. It is likely that this is the 
smallest ship that would carry an FEL. Eurthermore, and any effect calculated for this 
hull form would be reduced on any larger vessel. As such, if the system can be 
maintained in alignment on the DDG, it should be similarly possible on any larger ship. 
The worst-case sea state in which an engagement is likely to occur is sea state 6 (SS-6), 
from Table 3. Seas of 4 to 6m will be used in determining the steady state driving 
function for the model. 


Table 3: Table of Sea States. From [10] 


Sea 

State 

Description 

Wave 
Height (m) 

0 

Calm (glassy) 

0 

1 

Calm (rippled) 

0-0.1 

2 

Smooth (wavelets) 

0.1-0.5 

3 

Slight 

0.5-1.25 

4 

Moderate 

1.25-2.5 

5 

Rough 

2.5-4 

6 

Very rough 

4-6 

7 

High 

6-9 

8 

Very high 

9-14 

9 

Phenomenal 

over 14 
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Data from an underwater (UW) explosion was recorded during the shock trials of 
USS Winston S. Churchill and can be seen below in Figure 16. This will be 
approximated as a 45g, 4ms half sine pulse and applied to the model to explore how the 
system responds to shock loading from an underwater explosion. 
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Figure 16: Shock Trials Data for USS Winston S. Churchill. After [11] 


The DDG-51 hull form was entered into SFIIPMO [12], a program that tracks a ship’s 6 
degrees-of-freedom at various speeds and relative courses. MATLAB [13] was used for 
post processing of this data and polar plots of the resulting peak vertical acceleration and 
period for sea state 6 are below in 

Figure 17. In the figure, relative bearing to the seas is shown azimuthally, with 
180° indicating that the ship’s heading is directly into the seas, and the ship’s speed 
increasing radially from zero at the center to 30 knots moving outward. Considering the 
vibration isolation strategy that will be used, the worst inputs would be the highest 
amplitude for acceleration and the lowest period for that acceleration. As it turns out, 
these two conditions coincide at a single point, 30 knots with the ship’s heading directly 
into the seas. A worst-case steady state sinusoidal driving function of 0.12g in amplitude 
and a period of 5 seconds were determined. 
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* Relative Bearing of 180 is ships head into the seas 


Figure 17: Vertical Acceleration and Period for SS-6 Plotted Azimuthally vs. Ship’s Heading 

and Radially vs. Ship’s Speed 
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D. 


STATE-SPACE MODEL AND ANALYSIS 


MATLAB is used to perform a series of ealculations. First, it is used to eonvert 
the eigenvectors and values into a state-space model. The state-space model is a 
frequency domain model that is used to calculate transfer functions for the various 
components being tracked. An example transfer function is below in Figure 18, showing 
how the undulator and mirrors respond to input frequencies between 0.1 and 100 Hz. For 
illustrative purposes some of the various vibrational inputs to the system have been 
indicated in their relative frequency range. 



Wave Action 


Propulsion 

Auxiliaries 


10 '^.... . 

io’ 10® io’ 10^ 


frequency. Hz 


Figure 18: Example Transfer Function for Resonator Components 

The state space model is then converted into a time domain system model to allow 
for the computation of actual displacements and time histories in response to a driving 
function. An example motion plot is below in Figure 19; the figure shows the response in 
the vertical, (Z) direction, of a preliminary model that has four springs at the outer 
comers, from a continuous sine wave input. This time domain model also allows for the 
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comparison of motion between any of the components in the PEL system in any direetion 
and will be the primary design tool for optimizing various parameters to maximize 
system performanee. 


FEL mcKlet uz spring displacement vs time for 0.12g. 0.2 hz sinuslidal excitation 



Figure 19: Example Motion Plot 
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IV. SYSTEM DEVELOPMENT 


A. INPUT FREQUENCIES 

The typical spectrum of frequencies driving the FEL system onboard a ship is in 
Table 4. This set of frequencies lends itself naturally to an isolation frequency of 1 Hz 
because it is roughly a factor of 4 away from both of the largest amplitude sources of 
vibration, the waves and the propellers. 


Table 4: Environmental Vibrations 


Source 

Frequency Range (Hz) 

Wave Action 

<0.2 

Propulsion 

4-33 [14] 

Electrical Power 

60 

Auxiliary Machinery 

>60 


B, ISOLATION STRATEGY 

In order to attain the 1 Hz isolation frequency for the system, the spring constant 
must be calculated. The natural frequency of the isolation is f = where k is 

the spring constant in N/m, m is the suspended mass in kg, and f is the natural frequency 
in Hz; therefore, k = Tn^fm. Since the spring constant depends on the suspended mass, 
the baseplate design will have to be finalized prior to selecting it. 

Assume for the time being that the wave action frequency is ~0.1 Hz. If the 
isolation frequency is a decade away at 1 Hz, then there is sufficient separation that little 
to no energy transfer will occur between the excitation and the isolation. The same 
principle applies between the springs that are suspending the system and the system itself. 
With an -'1 Hz natural frequency for the springs, designing the baseplate for a -'10 Hz 
lowest distortion mode will place a decade of separation between it and the springs as 
well and minimize energy transfer between them. 
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c. 


BASEPLATE EVOLUTION 


1. Solid Plate 

The solid model was ereated with a ~160 ton, 10 eentimeter thick solid steel plate 
as its base. There was no particular justification for this dimension other than it seemed 
to be a reasonable starting point. The plate was then transferred to ANSYS and analyzed. 

It was determined that a 10 cm solid was inadequate to attain the 10 Hz or greater 
lowest distortion frequency desired. This simulation also revealed that the first distortion 
mode for this structure would be a bending mode and the second mode would be 
torsional; these mode shapes are shown in Figure 20. Several iterations later, it was 
determined that in order to achieve ~10 Hz, the baseplate would need to be 1 meter thick 
and weigh nearly 1600 tons. 



Figure 20: First Bending Mode (Left), First Torsional Mode (Right) 


2, Simulated I-beams 

A weight of 1600 tons is quite large for a single piece of equipment, so a better 
strategy for baseplate design would be necessary. A hand calculation revealed an I-beam 
array structure that has -95% of the bending stiffness of the solid plate, and has only 
-20% of the mass. Adding this into the ANSYS model as a lower density for steel and 
optimizing for first bending and torsional mode frequencies lead to an -150 ton structure 
that would perform as desired. There are, however, a couple of issues with this structure. 
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First, 150 tons is still heavier than it likely needs to be. Second, there is a reduction of 
torsional stiffness that arises from that use of I-beams that is not accounted for when 
simulating in this manner. 

3, Box Baseplate 

In continuing the development of the baseplate, a box structure was devised so 
that an accurate accounting of the internal structure could be achieved. The internal 
structure can be seen in Figure 21. To reduce overall system mass, the construction 
material was changed to aluminum. The model of this structure was set up so that a 
design study varying multiple parameters could be performed and an optimum could be 
selected. The thickness of the ribs, the thickness of the top and bottom plates, and the 
overall height of the structure could be varied individually while tracking the first 
bending and torsional mode frequencies and the mass of the system. 
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A parametric study was performed varying the rib, top and bottom plate 
thicknesses from 1 to 10 cm in 1 cm increments. Also, the box height was varied from 50 
to 150 cm in 10 cm increments. The results are in Figure 22, showing two points for each 
combination: one for the first bending mode and first torsional mode frequency versus 
the mass. From this plot, there emerges an easy selection for the system design, the 
rightmost of the three lowest points in the bending category. This case has a mass of ~55 
tons, a top and bottom thickness of 4 cm and a height of 1 m. 
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Figure 22: Box Baseplate Parametric Study Results 


This more practical structure led the design into another situation with problems. 
The intent was to minimize the system mass while attaining at least 10 Hz for the first 
bending mode, and as much separation between bending and torsional frequencies as 
possible. The heavier components, such as the linear accelerator cryo-modules and the 
beam dump, pull on the top plate at the mounts hard enough to deform it. This 
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deformation leads to excessive motion in the beam line and larger displacements for all of 
the components mounted to the bedplate. Figure 23 shows an exaggerated example of 
this type of deformation. 
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Figure 23: Distorted Surface from lusufficieut Top Plate Thickness 


To avoid significant top-plate deformation, the top plate needed to be at least 5 
cm thick. The extra thickness added enough additional stiffness so that the overall height 
of the baseplate could be reduced in order to lower the mass while maintaining sufficient 
performance. With the thicker top and bottom, 20 cm could be taken out of the height of 
the baseplate while keeping the 2 cm ribs and still maintaining the system performance. 
This resulted in the first bending mode occurring at 11.33 Flz, with the first torsional 
mode occurring at 21.12 Hz, and a mass of 58.2 metric tons. This combination is only 3 
tons heavier than the previous solution. 

D, SIMULATED SHIP STRUCTURE 

The starting point for the ship structure is shown in Figure 24. For this project, it 
is being assumed that the deformation of the ship’s hull will have little effect on a 
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resiliently isolated structure. Because of this, the simulated ship structure will be set up 
to minimize energy dissipation in it and ensure as much of the input as possible is 
transferred to the isolation system. This structure is constructed of beams that have 
effectively no mass and infinite stiffness so that it will not deform. In the center of the 
structure there is a mass node, an element with no volume that can be set to any mass 
needed. In this model it is set to 1000 times the mass of the rest of the system to ensure 
the motion of the suspended system does not cause a resonance in the mass node. This 
mass and beam combination will facilitate using the ship’s structure as a shaker table by 
applying motion to the mass node and translating that motion through rigid beams to the 
isolations. Four springs were used as a starting point in the vertical (Z) axis positioned 
-'1 m in from the outside edge of the baseplate, directly under the beam line, on the node 
line. The node line is a line that does not move when the system is excited in the first 
bending mode. In this case it is located '-'25% of the way in from both ends of the long 
(X) axis. Two more springs were used equidistant from the geometric center normal to 
the X axis, and one additional spring was used normal to the short (Y) axis. The 
idealized spring constants were calculated using the system mass determined earlier, the 
number of springs for a given axis, and the 1 Hz isolation frequency. 
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Figure 24: Initial Simulated Ship Structure 


Figure 25 shows the frequency response behavior of the mass node. This proves 
that the ship structure is sufficiently rigid, and the mass node is heavy enough to evenly 
convey the forcing function to the isolators without distorting or resonating it. If there 
were any nonlinearaties in this plot, the stiffness of the rigid beams or the mass of the 
mass node would have to be increased. 
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Figure 25: Mass Node Z Displacement vs. Frequency 


In this design, there is a significant difference in mass between the cryo module 
side of the PEL system and the optical resonator side. This difference generated a 
significant roll motion about the X axis which will need compensation. A first attempt to 
reduce this roll motion was made by moving the mounting location of the springs from 
below the beam line to the outer edge of the baseplate. Figure 26 shows the X direction 
motion of the left and right side mirrors relative to the undulator. From the figure it can 
be seen that the displacements have gotten about 10% smaller. This 10% reduction 
follows through most of the parameters being tracked by the model. 
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^ 0 '^ X delta displacement vs time for 0.12g. 0.2 hz sine input 



1.5- 


^ 0 '^ X delta displacement vs time for 0.12g. 0 2 hz sine input 



time, sec 


Figure 26: X-Direction Mirror Displacement Before (Top) and After (Bottom) Moving Mounts 


Consulting with the engineers at CSA-MOOG Engineering led to the selection of 
an off-the-shelf Firestone air spring to use as an isolator. These springs can support a 
static load of about 30 tons and have a spring constant of about 263 kN/m. Based on the 
mass determined for the baseplate, 1 Hz isolation requires a spring constant of 2.33x10^ 
kN/m. This means that 9 of these springs would be the closest match to the ideal case. 
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Considering that 9 springs is an odd number of real springs required, the model 
which only contained 4 springs was run with different spring constants on either side to 
simulate 9 springs and determine if the asymmetry would be beneficial or if 10 springs 
should be used. This step is being performed because it is far simpler to adjust a spring 
constant than to set up the ship structure to support more springs. This model was run 
with the greater spring constant on either of the long sides of the model. The roll motion, 
and the movement of the components attached to the baseplate, was worse with more 
springs under the lighter side. As expected, the greater spring constant under the heavier 
side further improved the system performance. The relative motion of the mirrors in the 
X-direction is shown in Figure 27, and shows a further reduction of relative motion 
between the two mirrors (forming the optical resonator) and the undulator. 
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Figure 27: X-Direction Mirror Displacement of Final Configuration 


Due to the success of the simulation of 9 springs, another model ship hull was 
created. This model had the 9 springs evenly spaced at the outer edge of the baseplate 
with 4 under the optical resonator and 5 under the heavier side. Since this model was 
only being excited in the vertical direction, the horizontal springs could continue to be 
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modeled as the ideal springs from the initial model. Figure 28 is the final ship structure 
on which the FEL system is mounted. It shows the vertical springs highlighted in green, 
the horizontal ones highlighted in red, and that more beams have been added to the 
structure to ensure deflections of the spring mounts are not introduced by asymmetries in 
spring location. 
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Figure 28: Final Simulated Ship Structure 
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V. SIMULATION RESULTS 


A, BASEPLATE MOTION 

The final configuration of the model was subject to both wave motion and shock. 
The resulting overall baseplate motion due to the sea input is shown in Figure 29 and 
shock is shown in Figure 30. The sea state 6 (SS-6) response resulted in about 60 mm of 
total travel with a period of ~5 s and some higher frequency transients that have periods 
of ~1 s. This travel approximately matches the linear response region of the selected air 
spring, and reinforces the Itkeliness that an actual system would perform as predicted by 
a linear model. The underwater (UW) explosion results are favorable as well. The 
springs are capable of a total motion of about 38 cm and under the worst-case input the 
system travels ~34 cm. The shock response produces motion with a period of ~1 s with 
no noticeable higher frequency components, and damps out at about 10% per second. 
This is desirable because there will not be an additional shock input as the spring reaches 
the end of its allowable travel and stops suddenly. 
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B. MIRROR ANGLE 

The angle of each mirror with respect to the central axis of the undulator was 
analyzed. Figure 31 shows an illustration of the mirror angle 0 that is being referred to in 
this section. The magnitude of mirror angular motion is summarized in Table 5. The 
following figures show the total angle between the undulator axis and the mirror axis. 
This angle is generated from rotations about the two non-symmetrical axes of the 
undulator. Figure 32 shows how the right mirror responds to the seas. It has the ~2.5 s 
period motion that is common in the steady state SS-6 responses in this model; it also has 
-0.5 s period components that are caused by the out of phase responses from the two 
source rotations. Figure 33 shows how the left mirror responds to the seas. It has the 
same -2.5 s motion as the right mirror but the higher frequencies are not as significant 
because the magnitude is about an order of magnitude greater. Figure 34 shows how the 
right mirror responds to the shock. It has the -0.5 s motion as well but it shows how this 
mirror takes -6 s to build up energy and begin to decay away. Figure 35 shows how the 
left mirror responds to the shock. It has the -0.5 s motion again but in this case there is a 
beating phenomenon that is likely due to the beam dump base structure not being 
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sufficiently rigid. The imbalances in the PEL system begin to become apparent in these 
plots. The left mirror, which is on the same side of the system as the beam dump, moves 
about one order of magnitude more than its twin on the opposite side. The table shows 
the peak values for the response to shock input and the steady state values for the 
continuous input once the startup transient subsides. The tolerance on the mirror angle is 
0.5 grad. Additional controls will be needed to control the mirror angle. 



Figure 31: Illustration of Mirror Angle 


Table 5:Mirror Angle Results 


Mirror 

SS-6 

UW explosion 

Left 

82 grad 

700 grad 

Right 

3.8 grad 

160 grad 
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Right mirror rotation 0.12g. 0.2 hz sine input 



20 30 

time, sec 


Figure 32: 


Right Mirror Angular Deflection in Response to SS-6 
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Figure 34: Right Mirror Angular Deflection in Response to UW Explosion 


Left mirror rotation 45g. 0.004 sec half-sine pulse 



Figure 35: Left Mirror Angular Deflection in Response to UW Explosion 
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c. 


MIRROR SHIFT 


The shift of each mirror in the plane normal to the central axis of the undulator 
was analyzed. An illustration of this shift (L) is shown in Figure 36; its magnitude is 
summarized in Table 6. Figure 37 shows how both of the mirrors respond to the seas. 
The shift is only a function of a single dimension and as such the short period motions 
occur at ~1 s vice the 0.5 s motion from the mirror angle results. Figure 38 shows how 
both of the mirrors respond the shock. The shock results show the same ~1 sec motion as 
the sea results and the beating caused by the beam dump are present again. The 
imbalances between the left and right sides of the FEL system continue to pose problems 
by causing the left side to move more than the right. The differences in the translation 
motions are not as severe as the angular ones; in this case the difference is closer to a 
factor of ~2. The tolerance on mirror shift is 50 pm; so that additional controls will be 
needed to control the mirror position. An active control system will have to reduce the 
effect of the shock by moving the mirrors by several millimeters in a about a tenth of a 
second to keep the system in operation. 



Figure 36: Illustration of Mirror Shift 
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Table 6: _ Mirror Shift Summary 


Mirror 

SS-6 

UW explosion 

Left 

1.2 mm 

7.8 mm 

Right 

0.8 

4.4 mm 


Mirror shift vs time 0.12g. 0.2 hz sine input 


1.5 



time, sec 


Figure 37: Mirror Shift in Response to SS-6 
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Figure 38: Mirror Shift in Response to UW Explosion 


D. OPTICAL PATH ANGLE 

The angle 0 between the line conneeting the two mirrors and the central axis of 
the undulator was also analyzed. Figure 39 shows an illustration of the angle 0 that is 
being referred to in this section. The seas deflect this angle 26 grad while the UW 
explosion deflects it 460 grad. Figure 40 shows the sea response of the optical path 
angle. It has the usual ~ 2.5 s motion as well as the ~1 s shorter period motion and a small 
amount of the -0.5 s motion in the first 5 s. The shorter period motion is not very 
significant in these results because of the larger amplitudes of the rotation. Figure 41 
shows the shock response of the optical path angle. It again displays the -3.5 s beating 
that is likely due to the beam dump’s weight and structure. SS-6 moves the optical path 
by -26 grad off axis, while the UW explosion moves it 460 grad. The tolerance on this 
motion is 1 mrad; the control system on the model is sufficient to maintain this parameter 
in tolerance. 
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Figure 39: Illustration of Optical Path angle 


Optical path angle vs time for 0 12g. 0.2 hz sinusoidal excitation 



Figure 40: Optical Path Angular Response to SS-6 
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Optical path angle vs time for 0.12g. 0 2 hz sinusoidal excitation 
500 



Figure 41: Optical Path Angular Response to UW Explosion 


E. OPTICAL PATH SHIFT 

The motion in the plane normal to the central axis of the undulator of the line 
connecting the two mirrors was analyzed. Figure 42 shows an illustration of the shift L 
that is being referred to in this section. Figure 43 shows the sea response of the optical 
path shift. It is only slightly affected by the ~1 s motion before it transitions to the steady 
state ~2.5 s period motion. Figure 44 shows the shock response of the optical path shift. 
It is done with any transient behavior in the first 2 s and then the amplitude decays away 
at -10% per second. The optical path is shifted -0.4 mm off axis by SS-6 and -2.3 mm 
by the shock. The tolerance here is 0.6 mm; the system will be out of tolerance for about 
-15 sec during an UW explosion without active controls, but the passive controls in this 
model will be sufficient during operations up to SS-6. An active control system will have 
to reduce the effect of the shock by moving the mirrors a few millimeters in about a tenth 
of a second to keep the system in operation. 
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Figure 42: Illustration of Optical Path Shift 
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Optical path shift vs time for 45g. 0.004 sec half-sine pulse 



Figure 44: Optical Path Shift Response to UW Explosion 


F. RESONATOR LENGTH 

The change of length of the resonator was analyzed. Figure 45 shows an 
illustration of the length L that is being referred to in this section. Figure 43 shows the 
resonator length behavior in SS-6 waves. The 1 s motion persists through the first 10 s of 
operation and then the steady state behavior is established. Figure 44 shows the resonator 
length behavior during an UW explosion. The high frequency behavior subsides in the 
first 2 s and then the ~2.5 s motion decays away at -10% per second. Under SS-6 input 
the resonator grows and shrinks by about 90 pm. Under UW explosion the resonator 
grows and shrinks by about 500 pm but returns to the SS-6 value in about 30 s. The 
tolerance on this motion is 0.5 pm and is essential for maintaining the FEL in operation. 
Further control will be needed but this model provides an operational envelope in terms 
of period and amplitude in which the active control system must operate. 
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Figure 45: Illustration of Resonator Length 


Change in resonator length vs time for 0 12g. 0 2 hz sinusoidal excitation 



Figure 46: Resonator Length Response to SS-6 
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VI. CONCLUSION 


A, CONCLUSION 

The tolerances and worst-case result for sea state 6 (SS-6) waves and under water 
(UW) explosion are summarized in Table 7, with out of tolerance values in red. Three 
out of the 5 figures of merit developed in this project are not in tolerance for the worst- 
case scenarios chosen. The other sea-states were also explored to show how the system 
behavior changed with weakening seas as shown in Table 8. All of the parameters being 
monitored show significant decrease with sea-state but only mirror shift changes to being 
in tolerance and that only happens in sea-state 3. 


Table 7: Tolerances and Simnlation Resnlts 


Dimension to be Controlled 

Tolerance 

Model 

Prediction for 
Sea State 6 

Model Prediction 
for Under Water 
Explosion 

Mirror Angle 

±0.5 grad 

80 grad 

700 grad 

Mirror Shift 

±50 gm 

1.4 mm 

8 mm 

Optical Path Angle 

±1 mrad 

26 grad 

460 grad 

Optical Path Shift 

±0.6 mm 

400 gm 

2.3 mm 

Mirror Separation 

±0.5 gm 

90 gm 

500 gm 


Table 8: Simnlation Resnlts for Varied Sea States 


Sea State 

Dimension ^- 

6 

5 

4 

3 

Mirror Angle 

80 grad 

34 grad 

2.2 grad 

3.4 grad 

Mirror Shift 

1.4 mm 

480 gm 

90 gm 

48 gm 

Optical Path Angle 

26 grad 

12 grad 

2 grad 

1.1 grad 

Optical Path Shift 

400 gm 

160 gm 

33 gm 

16 gm 

Mirror Separation 

90 gm 

38 gm 

7 gm 

3.8 gm 
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The passive control system laid out in this project is not able to maintain the PEL 
system components in alignment sufficiently to allow for its operation. An active control 
system will be required. This system does, however, reduce the amplitude of motions 
that the active control system must be able to control. 

Typical active controls systems operate in the kHz range and all of the motions in 
this system that need to be controlled are in the 0.5 to 5 Hz range. This, combined with 
the reduced amplitudes of motion from the passive isolation system, make it likely that 
the FEE system can be operated successfully onboard ships. 

The goal of this project was not to develop a passive isolation system that would 
enable FEE operation on board a ship. It was to develop a system that would reduce the 
amplitude of the vibrations left to active controls to a reasonable level and to establish the 
operating envelope for that system. Both goals have been successfully achieved. 

B. FUTURE WORK 

To improve this model, several steps can be taken in the future. Adding active 
controls to the air springs, such as a voice coil actuator, would greatly improve 
performance and would be closer to the system that will be installed on a ship. 
Furthermore, a more complete model of the ship’s hull could be added that would allow 
for the addition of the effect of hull bending on the model. Also, there are many FEE 
system components such as RE power and the cryo plant that were assumed not to be on 
the baseplate. Either adding them onto the baseplate or modeling their connections will 
be important in a better model. Finally, once the shipboard FEE system design matures 
to the point that no further major structural changes are going to be made, that model 
should be incorporated into this project. 
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